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Abstract 

The Lorentz-invariant nuclear lagrangian of Furnstahl, Serot and Tang (FST) is discussed. The 
FST lagrangian is derived in terms of an effective field theory and exhibits a nonlinear realization 
of chiral symmetry SU (2)^ x SU (2)#. The relevant degrees of freedom are nucleons, pions and the 
low-lying non-Goldstone bosons: isoscalar scalar (<r) and vector (w) mesons, and isovector vector (p) 
mesons. The terms in the lagrangian are organized by applying Georgi's naive dimensional analysis 
and naturalness condition. As a consequence all coupling constants in theory are dimensionless 
and of order unity. 
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I. INTRODUCTION 



The effective field theory (EFT) technique allows to construct in a controlled manner, 
below a characteristic energy scale, the most general lagrangian consistent with relevant 
degrees of freedom and symmetries of an underlying theory. In nuclear physics, the EFT 
method relies on the symmetries of QCD to construct the effective lagrangian. The main 
component of this construction is the chiral SU(2)l x SU(2)r symmetry. This almost 
perfect symmetry is spontaneously broken to its vectorial subgroup SU(2)v=l+r with the 
appearance of pseudo-Goldstone bosons (pions). 

One of the recent attempts at formulating EFT for finite nuclei and nuclear matter is 



the generalization of Wa 



ecka quantum hadrodynamics (QHD) proposed by Furnstahl, 
Serot and Tang (FST) The FST lagrangian is derived by expansion in powers of 

the lowest lying hadronic fields and their derivatives. The relevant degrees of freedom 
are nucleons, pions, isoscalar- vector fields (tu-mesons) , isoscalar-scalar fields (a-mesons), 
and isovector- vector fields (p- mesons). In this lagrangian the chiral symmetry is realized 
nonlinearly using a standard procedure of Weinberg [5J], and Callen, Coleman, Wess, Zumino 

n nn 

|6|, (WCCWZ). The terms of the la grangian are organized by applying Georgi's [7|, |8[ naive 
dimensional analysis (NDA) and a "naturalness" condition. 

The framework of EFTs and Georgi's naive dimensional analysis are detailed in Sect. [Til 
In Sect. II I II the nonlinear realization of chiral symmetry and the FST chirally symmetric 
effective lagrangian are discussed. The Dirac-Hartree approximation of the FST lagrangian, 
i.e., a treatment of the lagrangian at the level of classical meson fields and valence nucleons 
is shown in Sect. IIVI A short summary is given in the last section. 

II. EFFECTIVE FIELD THEORIES 



A. Principles of the EFTs 
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An essential idea underlying the effective field theories (EFTs), see e.g. Ref. 
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15], is relevant to the appearance of disparate characteristic energy scales, 
E <C E , in quantum field theories. Suppose that we are interested in physics at lower scale 
E, then we can choose a cut-off scale A at or slightly below Eq and divide the generic fields 
into two parts: a low- 4>l and a high-energy (pu (<P = 4>l + 4>h), according their momenta 



are smaller or greater than A 



<j> L {k): \k\<K, (j> H {k): \k\>K. (1) 

The effective lagrangian is obtained by path integrating over the high-energy part <pn m the 
generating functional Z 

Z = J [d(j) L ] [d<p H ] j$^t<l>M = J [d(j) L \ e iJ^B(M t (2) 

where 

J d 4 x£ efr (0 L ) = -ilnj [dc/> H ] e*f^^M. (3) 

This defines the procedure of eliminating the high-energy degrees of freedom referred 
to as "decimation". The next step is to write £ e fr in terms of local operators O^^l), 

oo 

£eff(0r) =X>(A)O i (&0, ( 4 ) 

i 

where ft (A) are the coupling constants absorbing the contribution of the integrated out 
high-energy degrees of freedom 

In Eq. (j3J) the effective lagrangian is represented by a infinity series of interactions that 
involve the relevant degrees of freedom and satisfy the assumed symmetries of the underlying 
high-energy theory. In order to make this procedure useful we need some dimensional 
analysis. In units h = c=l the action S e ff = ^d^xC^ (4>l) is dimensionless. If an operator 
Oi has dimension Si, [Oj\ = [m] 6 ' = 6i, then gi has dimension [ft] = 4 — 8\ and we can define 
dimensionless coefficients q = A^ _4 ft, which additionally are assumed to be "natural", i.e. 
of order 0(1). For a process at scale E, we can estimate dimensionally the magnitude of 
the i'th operator in the action as 

J d A xOi ~ E 5 '-\ (5) 

so that the i'th term is of order 

Now we can see that at energies below A, the behaviour of the different operators is deter- 
mined by their dimension. If 8{ < 4 the operator is more and more important when E — > 0, 
and is termed relevant. Similarly, if Si > 4 the operator is less and less important, and is 



termed irrelevant. An operator with Si = 4 is equally important at all scales and is called 
marginal. 

At energies much below A, corrections due to the irrelevant (non-renormalizable) parts 
are suppressed by powers of E/A and the effective lagrangian is able to describe the low- 
energy physics. The accurate procedure which connects the order of the expansion in powers 
E /A with the terms in the effective lagrangian that need to be included at that order is called 
"power counting". 



B. Naive dimensional analysis 

There are at least two relevant energy scales in nuclear physics: the pion-decay constant 
f n ~ 93 MeV and the larger scale A ~ 4:irf n ~ f GeV, which characterizes the mass scale 
of physics beyond Goldstone bosons. Using a naive dimensional analysis (NDA) proposed 
by Georgi and Manohar EJ, for assigning the LECs of appropriate sizes, the effective 
lagrangian describing interactions of nucleons N(x), pions tt(x), and non-Goldstone bosons 
(scalar (f)(x) and/or vector V(x) mesons) takes a general form 



DC 



= £^ (A) ' (8) 

A=0 

where T is a product of Dirac matrices, T> a covariant derivative, a pion mass (treated 
as derivative) and c n d p b the dimensionless LECs which are assumed to be natural, of 0(1). 

In Eq. (JHJ) the interactions are grouped in sets £' A ' of common index A = & + d + b — 2, 
according to ((71), each of them carries a factor of the form (1/A) A . This formula has profound 
implications if we invoke chiral symmetry, Ref. If]]. For strong interactions (in absence of 
external gauge fields, e.g. photons) chiral constraint guarantees that A > and the large 
scale A~47r/ 7r ~l GeV does not occur with positive powers in Eq. (J7J). 
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III. EFFECTIVE CHIRAL LAGRANGIAN FOR NUCLEI 



A. Nonlinear realization of chiral symmetry 

Let us briefly collect the basic ingredients for considering the chiral effective lagrangian. 
In the chiral limit, where Nf = 2 or 3 of quarks are massless (u, d and possibly s), the 
underlying QCD lagrangian is invariant under a global group 

U(N f ) L x U(N f ) R ~ SU(N f ) L x SU(N f ) R x U(l) v x £7(1)^. (9) 

However, at the quantum level, due to the axial anomaly, the U(1)a symmetry is broken 
that, e.g., leads to nonzero mass of 77' (958) meson even in the chiral limit. In the hadronic 
world, the chiral group G = SU(Nf) L x SU(Nf) R is spontaneously broken in the vacuum 

SU(N f ) L x SU(N f ) R x U(l) v — > SU(N f ) v=L+R x U(l) v , (10) 

to the vectorial subgroup H = SU(Nf)y, either of isospin when Nf = 2 or flavor SU(3) 
when Nf = 3. The preserve vector group U(l)y is realized as baryon number conservation. 
According to Goldstone's theorem the number of Goldstone fields is the dimension of coset 
space G/H, which is the number of generators of G that are not also generators of H. In 
our case of chiral symmetry dim G/H = Nj — 1 and we can identify Nj — 1 pseudoscalar 
Goldstone bosons, ip = 7r[, K, 77], with the pions for Nf = 2, plus the kaons and an eta meson 
for N f = 3. 

The nonlinear realization of spontaneously broken chiral symmetry, denoted by WCCWZ, 
was suggested by Weinberg f5] and developed further by Callan, Coleman, Wess, and Zumino 
jfj]. In the WCCWZ formalism the Goldstone bosons (p, being coordinates of the coset space 
G/H, are naturally represented by elements £(x) = £(<p(x)) of this coset space. The chiral 
symmetry is defined by specifying the action of G on the representative £(<£>), with a canonical 
choice of coset representative this transformation takes the form 

£(</?) 9R^)h\g, if) = h{g, <p)({<p)gl (11) 

where g = (gL,g R ) G G. The equality in (TTT!) is due to parity and it defines the so-called 
compensator (field) h(g,<p(x)) G H. Its dependence on the Goldstone boson fields ip(x) is a 
characteristic feature of the nonlinear realization of chiral symmetry. 
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Let us restrict ourselves to the two flavor case Nj = 2 with the isotriplet of pions collected 
in a 2 x 2 special unitary matrix 



§t*=§ , (12) 

yV2n- -7T° J 

where f are Pauli matrices. Applying an exponential parametrization a coset representative 
£(x) = ^(tt(x)) can be written as 

U(x) = e(x) = exp {2m{x)/U) , (13) 

where f w ~ 93 MeV is the pion-decay constant. The isospinor nucleon field is represented 
by a column matrix 

\n{x) J 

where p(x) and n(x) are the proton and neutron fields, respectively. The relevant low-lying 
non-Goldstone bosons are an isovector- vector p(770) meson p^(x) = p^x) ■ |r, an isoscalar- 
vector meson u;(782) represented by a vector field V^x) and an effective isoscalar-scalar field 
S(x) (a meson) is described by the shifted field (Kx) = S — S(x), where S is the vacuum 



expectation value of the scalar field S, see Ref. [17| . The to meson is needed to describe the 
short-range repulsion and effective field is included to incorporate the mid-range attraction 
of the NN interaction. 

Following WCCWZ, the nonlinear realization of the chiral symmetry for the mentioned 
above degrees of freedom are 

9r€(x) h ] {g,-K(x)) = h(g,n(x)) £(x) g{, (15) 
N(x) h(g,7r(x))N(x), N(x)-^N(x)h\g,Tr(x)), (16) 
Pr(?) h(g,7t(x)) p^x) rf(g,ir(x)), (17) 

where g = (gL,gn) £ SU(2)l x SU(2)r. The matrix field U(x) transforms linearly under 
chiral transformations 

U(x) ^ g L U(x)gl (18) 

and the other degrees of freedom, neutral-isoscalar fields V^(x) and <f>(x), can be treated as 
the chiral singlets. 
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We can also define an axial vector field a M and a polar vector field 

a , = ~(e t s / .e-^e t ) = 4' ( 19 ) 

v, = - l -^ d ^ + ^) = vl, (20) 

where the hermiticity follows from c? M (£'£) = = <9 M (££^). Under the chiral symmetry the 
transformation of is homogeneous 

/ic^/i 1 ", (21) 

whereas that of is inhomogeneous 

v„ h^-id^rf. (22) 

In fact is the connection on the coset space and with it we can construct the covariant 
derivatives on this space. For example, since transformations ( 1T6]) and ( fTTj) are not only 
nonlinear but also local, via compensator field h(g,n(x)), that requires the introduction of 
the chirally covariant derivatives for nucleon and rho meson fields: 

Dfj,N = (d^ + iv^N, D^p v = d^p u + i[v„, p v ], (23) 

which transform covariantly under the chiral group. Also, the curvature (strength) tensor, 
v^v, associated with the connection can be expressed in terms of axial vector fields, a^, as 

V = dpV,, - d v v^ + i[Vft, v v ] = -i[a,t, a v \. (24) 

The covariant derivative of rho meson can be used to construct the covariant field tensor 

PHV = Vy.Pv ~ VvPy, + ig P [Pn: pu], (25) 

which is the last ingredient needed to build the FST effective chiral lagrangian. 



B. FST effective chiral lagrangian 



The 
Ref. 



18 



ow-energy effective lagrangian of Furnstahl, Serot and Tang (FST) [3j, |4j, see also 



19 



20], incorporates the symmetries of QCD: Lorentz invariance, parity invariance, 



nonlinear realization of chiral SU(2)l x SU(2)r symmetry, this lagrangian is also invariant 
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under the electromagnetic U(l) em and isospin SU(2) groups. The FST lagrangian is ex- 
panded in powers of fields and their derivatives in the procedure of power counting with 
index A = ^ + d + b, where n is the number of nucleon fields, d is the number of deriva- 
tives and b is the number of non-Goldstone boson fields in each term. Taking as the large 
energy scale A in Eq. (JTj) the nucleon mass M = 939 MeV, we may write the effective chiral 
lagrangian through quartic order (A < 4) as the sum 



C$(x)+£§(x)+£ 



(4) 

EM 



X) 



(26) 



The part of lagrangian involving nucleons takes the form 



N Y {idf, -v p - g p p p - g v V^) + # A 7 M 75<V — (M — g s <p) 

fp9p -N Puu a^N - hElNV a ^N - ^-Nv uv a» v N 
4M M M M 



N 



AM 



+ 



(27) 



where <r lu , = i[y li ,'y v \, V pv = d p V u — dyV^ is the covariant tensor of the ui meson, g A ^1.26 is 
the axial coupling consta nt, q n , f p , and g Y , f v are vector and so-called tensor couplings for 



p and uo mesons, see Ref. 21], g s is a Yukawa coupling for the effective scalar field 0, and 



is the coupling for higher-order irN interaction. The ellipsis represents redundant 



or tiny additional terms with irN and nir interactions, which have been omitted in the FST 
lagrangian. 

The mesonic part of the effective lagrangian is 

f 2 

J 7T 



\ (l + c*i^) d^dy + ^ tr (d,Ud^) 
~\ tr (P^ u ) ~ \ (l + «2^) V^W" - g, 2 % 



7V7T o 

nip 



tr {pu„ir) 



+ (l + r l P jf)m 2 p tr (p p p») - m 2 J 



(28) 



2 3! M 4! M 2 



where m v = 782 MeV, m p = 770 MeV, and m s are u, p and a mesons masses, g p7T7T is the 
coupling of pw interaction, which (assuming vector-meson dominance) is g pn7T = g p . The 
trace "tr" is in the 2x2 isospin space. 
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The electromagnetic interactions are described by 

4 4 M = -\F^F^-\eN^{i + T ^) NA ^- ^F, v N\a^N 
e — ~ — ~ ,2 



2M 2 



N llx {(3 s + p v r 3 )Nd v F^ - 2ef*A» tr (v M r 3 ) 

_,„. tr(r 3P n + |H+---, (29) 

where A p is the electromagnetic field, F pu is the electromagnetic field-strength tensor. Ac- 
cording to vector- meson dominance and phenomenology one can find that g 7 = 5.01. The 
lagrangian C^ M is invariant under the U(l) em group. The composite structure of the nucleon 
is included through an anomalous moment A = ^A p (l + t 3 ) + ^A n (l — r 3 ), with A p = 1.793 and 
A n = —1.913 the anomalous magnetic moments of the proton and the neutron, respectively. 
The ellipsis represents redundant terms of 0(e 2 ). 

The effective chiral lagrangian Eq. (1261) at a given order contains certain parameters 
that are not constrained by the symmetries, the so-called low-energy constants (LECs). 
Apart from the isoscalar (/3 B ), isovector (j3 v ) electromagnetic form factors and the tensor 
coupling for p meson (f p ), which are fixed from the free-space charge radii of the nucleon, 
the remaining thirteen LECs { ||, f|, f£, 771, r/ 2 , r) p , k 3 , k 4 , Co, , /v, "i, "2} have to be 
determined from experimental data. The LECs are defined applying the naive dimensional 
analysis so that they are assumed to be of order unity, i, e. "natural" . 



IV. DIRAC-HARTREE APPROXIMATION 



The mean-field approximation (ignores) dismisses all quantum fluctuations of the meson 
fields and treats them by their expectation values. Assuming the time reversal invariance 
the spatial components of the vector field vanish and we can define scaled mean meson fields 
(potentials) by including couplings: W(r) = g v V (r), $(r) = g s (p Q (r), R(r) = g p b (r) and 
A{r) = eA (r). Since the nuclear ground state has a well-defined charge, only the neutral 
rho meson field (denoted by b ) has been used, also since the ground state is assumed to have 
well-defined parity the pseudo-scalar pion field does not contribute in this approximation. 

If we restrict consideration to static nuclear systems the Dirac equation with eigenvalues 
E a and eigenf unctions ip a ( r ) is, see Ref. j^J, 

hip a (r) = E a tp a {r), J d 3 xiPl{r)ip a {r) = 1, (30) 
9 



with 



h(r) 



-ia • V + W(r) + \r^R{r) + (3 M - ${r) + \{l + r 3 )A(r) 



2M 

i 



pa- 



f±r 3 VR(r) + f v VW(r) 



2M 2 



(A + p v r 3 ) AA(r) 



\(5cx ■ VA(r) 



2M 

where f3 = 7 , a = 707. 

The mean-field equations for $(r), W(r), -R(r) and A(r) are 



(-A + m 2 ) #(r) 



^ 2 Ps(r)-^^ 2 (r) 



2ML 

Oil 



2M 



#(r) 
, 2 



K3 k 4 0Q) 
2 3! M 



™f W 2 (r) + 9 -^^R 2 (r) 



9v 



2 W 2M d 



«2« 



2M#. 



(-A + m 2 ,)^) = ^ 



(V#(r)) +2#(r)A$(r) 
-(W(r)) 2 , 

PB(r) + ^V.(fpI(r)) 



+ 



'/1 
M 
3# 7 



772 



2 M 

V<2>(r) • VW(r) + tf>(r)AW(r 

Pchg(r), 



(-A + mJ)J2(r) = i<? 2 L(0 + A V .(r^(r)) 



2M 

-^^^(r)-^ Pchg (r) 



- AA(r) = e 2 p chg (r). 

Assuming spherical symmetry and parity conservation the eigenfunctions of Dirac 
tion (I30I) (the positive-energy spinors) can be written as 



ip a (r) = ip nKm t(r) 



1 [G a {r)/r] $ K 



- [F a (r)/r] <P_ Kr 
dr(\G a (r)\ 2 + \F a (r)\ 2 )=l 



Ct, a = {n,K,t}, 
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where $ Km = J2m e m s (^ m el m s\j m ) Y £,rn e (^)Xm s are spin spherical harmonics, n is the prin- 
cipal quantum number, k is a nonzero integer uniquely determining j and £ through 
k = (2j + 1)(£ — j) and ( t is a two-component isospinor labeled by the isospin projec- 
tion t = \ for protons and t — — \ for neutrons. The radial equations for upper (G) and 
lower (F) components become 

d + ^) G ^ r ) - [ E « - U ^ r ) + U 2(r)]F a (r) - U 3 (r)G a (r) = 0, (38) 

ilTr ~ ^) Fa{r) + ^ " Ul{r) " U *W\ G *W + u *W F *<r) = °' ( 39 ) 

where single-particle potentials are given by 

Ul{r ) = W{r ) + tajR(r) + (f o + i)A(r) + ^2 (/3s + 2t a /3 v )AA(r), (40) 

C/ 2 (r) = M — <2>(r), (41) 
dW(r) , dR(r) 



[±(A P + A n ) + t a (X p - X n )] } . (42) 



The various densities that appear on the r.h.s. of the meson equations for spherically sym- 
metric systems are defined as follows: 

occ occ _ . 1 

Ps(r) = 5>«(rty a (r) = E ^!L+ (Gftr) - F a 2 (r)) , (43) 



a 

occ 



occ occ ~ . 1 

PbW = ^^(r)i)9a.f^ a (r) = E "i^T 2 G a (r)F a (r), (45) 

a a 
occ occ ~ . 1 



a a 
occ occ 



Pl(r) = X>£(r)irs/?a.nMr) = E 2 G a(r)F«(r), (47) 

a a 

where the summation superscript "occ" means that the sum runs only over occupied (va- 
lence) orbitals up to some value of n and k. The quantum numbers are denoted by 
{a} = {a; m} = {n, k, t; m}. 
The charge density is given by 

Pchg(r\) = p d (r) +(>m(r), (48) 
11 



where the "direct" nucleoli charge density is 



Pd(r) = P P (r) + ^V- [?p T a (r)} + ^ % Ap B (r) + (3 Y Ap 3 (r)} 



and the vector mesons contribution is 



Pm(r) 



-AR(r) 



2M 2 



1 



(49) 



-AW(r). 



Here the "point" proton density and nucleon tensor density are given by 

occ 

p p (r) = 1 1 J2^l(r)(l + T 3 )Mr) = l[PB(r)+p 3 (r)}, 

a 

occ 



(50) 



(51) 
(52) 



respectively, where A is the anomalous magnetic moment. Thus the spherical nuclear ground 
state with the presence of time reversal symmetry is described by coupled, one-dimensional 
differential equations that may be solved by an iterative procedure. Once the solution has 
been found, the total energy of the system is given by 



OCC p 

E = Y,E a {2j a + l)- / d 3 xU„ 



(53) 



where 



U„ 



ml $ 3 



+ * Af>A+ g* M -12 2 IV 



2M 
k 3 k 4 $ 



P3 + 4V.(^) 



AM gj 



<P 



AM 
AgjM 



<P i ml 

2 



1 



4!^ 



ywy 



(54) 



One of the most prominent observables, the binding energy of a system of A = Z + iV 
nucleons is defined by 

E B = E - Ecu - AM, (55) 

where -Ecm is the center-of-mass (cm.) correction which can be estimated nonrelativistically, 
e.g., its an empirical estimate given by Reinhard[22] is Ecm ~ 17.2A -0 ' 2 MeV . An older 
estimate from the harmonic oscillator shell model is Ecm ~ §41t4~ 1 / 3 MeV. 
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The mean-square radius of the charge distribution, with the (cm.) motion correction, is 
given by 

<Och g = (r 2 ) - (56) 



where 



(r 2 ) = ^J d 3 xr 2 p chg (r), (P 2 M ) = 2AME CM . (57) 

Since the additional nonrenormalizable interaction between the nucleon and electromagnetic 
field were included in £em, Eq. (129|) . the charge density p c h g automatically contains the 
effects of nucleon structure, and it is unnecessary to introduce an ad hoc form factor in 
formula (|56|) . 



V. SUMMARY 



One of the major challenges in nuclear physics is to establish a connection between nuclear 
dynamics and the fundamental QCD. The chiral effective field theories are considered to offer 
a natural and useful framework for this purpose. Thanks to the implementation of nonlinear 
realization of chiral symmetry, Georgi's naive dimensional analysis and the "naturalness" 
condition, the FST approach is the extension of Walecka's hadrodynamics and may be used 
in nuclear physics to cross the border from QCD to a nuclear theory. 
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